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An experimental study has been conducted to assess the effects of compression pad cavities on the aeroheating 
environment of the Project Orion Crew Exploration Vehicle heat shield. Testing was conducted in Mach 6 and 10 
perfect-gas wind tunnels to obtain heating measurements in and around the compression pads cavities using global 
phosphor thermography. Data were obtained over a wide range of Reynolds numbers that produced laminar, 
transitional, and turbulent flow within and downstream of the cavities. The effects of cavity dimensions on boundary- 
layer transition and heating augmentation levels were studied. Correlations were developed for transition onset and 
for the average cavity-heating augmentation. 
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Nomenclature 

cavity diameter, m 

cavity aspect ratio 

cavity depth, m 

total enthalpy, J /kg 

cold-wall (at 300 K) enthalpy, J /kg 

ratio of measured-to-Fay-Riddell-theory heating 

coefficients 

ratio of measured heating coefficients with and 
without cavity 

boundary-layer-edge Mach number 

freestream Mach number 

freestream pressure, Pa 

heat flux predicted using Fay-Riddell theory, 

W/cm 2 

maximum vehicle radius, m 

freestream unit Reynolds number, 1 /m or 1 /ft 

boundary-layer-edge momentum thickness 

Reynolds number 

nose (spherical cap) radius, m 

radius at aftbody shoulder, m 

radius at tangency point of spherical cap and 

shoulder, m 

freestream temperature, K 

freestream velocity, m/s 

distance along vehicle symmetry (roll) axis, 

m or in. 

distance in vehicle yaw-plane, m or in. 
distance in vehicle pitch-plane, m or in. 
angle of attack, deg 
boundary-layer thickness, m 
boundary -layer momentum thickness, m 
freestream density, kg/m 3 


I. Introduction 

T HE Crew Exploration Vehicle (CEV) concept was defined by 
NASA’s Exploration Systems Architecture Study [1]. This 
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study was conducted in 2005 to develop requirements for crew and 
cargo launch systems to support lunar and Mars exploration 
programs as well as for access to the International Space Station 
(ISS). Several versions of the Orion CEV are planned that will 
provide transportation first to the ISS, then to the moon, and finally to 
Mars. 

The crew module of the CEV (Fig. 1) has a configuration that is 
outwardly similar to the Apollo command module: a spherical- 
segment heat shield joined by a small toroidal section to a truncated- 
cone- shaped crew compartment. The Orion CEV, however, will be 
considerably larger than Apollo with a maximum heat-shield 
diameter of ~5 m (current configuration) vs 3.912 m for Apollo. 
This larger size will allow transport of up to six crew members on ISS 
missions or up to four crew members on lunar missions. 

In a proposed design similar to that of the Apollo program, the 
CEV crew module will be attached to the service module by six 
compression pad and tie-rod assemblies (Fig. 2) spaced at 60 deg 
increments around the heat shield of the vehicle. Because the 
compression pads and the heat shield will be made from different 
materials, they will ablate at different rates during reentry and so at 
different points along the reentry trajectory the compression pads 
will be either recessed below, or protruding above the heat shield. As 
evidenced from examination of recovered Apollo heat shield (Fig. 3) 
the compression pads, whether recessed or protruding, have the 
potential to cause large, albeit localized perturbations to heating 
levels as well as to promote boundary-layer transition. 

To obtain information to support the design of the CEV heat shield, 
a parametric investigation of compression pad dimension effects on 
heating has been performed on a generic CEV design through testing 
in the Langley Research Center 20 Inch Mach 6 and 3 1 Inch Mach 10 
Air Tunnels. In these tests, global heat-transfer data on the pads and 
heat shield were obtained using the Langley Research Center ther- 
mographic phosphor imaging technique. These tests supplement 
earlier efforts [2] in which models with the same cavity diameter but 
with larger cavity depths (hence smaller aspect ratios, D/H) were 
tested. Additional data were also drawn from [3] in which various 
cavity sizes and configurations (including tension ties and annular 
attachment points) were tested. Supporting information was obtained 
through computational fluid dynamics predictions for “smooth- 
body” (no-compression) pads performed at the wind-tunnel test 
conditions to provide boundary-layer properties for use in correlation 
of the cavity transition and heating data. 


II. Experimental and Computational Methods 

A. Compression Pad Wind-Tunnel Model Design 

Wind-tunnel models of the CEV vehicle with compression pad 
cavities were slip-cast in silica-ceramic and coated with a thermog- 
raphic phosphor compound [4], The model outer mold line (OML) 
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Fig. 1 NASA CEV (conceptual artwork). 



Fig. 2 Generic compression pad/tie-rod assembly for CEV. 


was a generic smooth-body representation based on early Project 
Orion design cycles (Fig. 4). The model scale was 0.03556 (based on 
an assumed 5.00 m full scale at the time of the model design) that 
resulted in baseline 7-in.-max-diam models, although an additional, 
larger 8-in.-diam model was built to obtain naturally turbulent data. 
Compression pads were located on each model at 60 deg increments 
around the circumference of the heat shield at a radial position of 
81% of the maximum radius (Fig. 5). The baseline pad diameter was 
0.285 in. on the model, which represented full-sized pads of 8.015 in. 
based on geometric scaling. Because the ablation/recession response 
of the pads and surrounding heat shield are not yet defined and the 



Fig. 3 Effects of compression pad heating on Apollo command module 
heat shield. 


compression pad design has not yet been finalized, models with 
cavity depths of 0.00625 in. (C0-D1 model), 0.0125 in. (Cl-Dl 
model), and 0.0250 in. (C2-D1 model) were fabricated. Cross- 
sectional views of these various cavity geometries are shown in 
Fig. 6. Also shown in this figure is a larger 0.0500 in. depth model 
(C3-D1 model) tested in the previous phase of the program [2]; note 
that the IX and 2X model designations assigned in that test have been 
redesignated as C2-D1 and C3-D1 for consistency with the current 
work. A final model (C1-D2) with a depth of 0.0125 in. but with half 
the nominal diameter (0.1425 in.) was also tested and is shown in this 
figure. 


B. Facility Descriptions and Test Parameters 

1. NASA Langley Research Center 31 Inch Mach 10 Air Tunnel 

The NASA Langley Research Center 3 1 Inch Mach 10 Air Tunnel 
(Fig. 7) is a conventional blowdown facility in which heated, dried, 
and filtered air is used as the test gas; a detailed description of this 
facility can be found in [5]. The tunnel has a two-dimensional, 
contoured nozzle that opens into a 3 1.0 by 3 1.0 in. test section. The 
tunnel is equipped with a side-mounted injection system that can 
transfer a model from the sheltered model box to the tunnel centerline 
in less than 1 s. Run times of several minutes are possible in this 
facility, although for aeroheating tests, run times of only a few 



Fig. 4 CEV wind-tunnel model dimensions. 
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C0-D1 model, D/H = 45.6 
Height = 0.00625 in. 
Diameter = 0.285 in. 


C1-D1 model, D/H = 22.8 
Height = 0.0125 in. 
Diameter = 0.285 in. 


C2-D1 model, D/H = 11.4 
Height = 0.0250 in. 
Diameter = 0.285 in. 


C3-D1 model, D/H = 5.7 
Height = 0.0500 in. 
Diameter = 0.285 in. 


C1-D2 model, D/H = 11.4 
Height = 0.0125 in. 

Diameter = 0.1425 in. 

1 i 

Fig. 6 Compression pad cavity cross sections. 

seconds are required. The nominal reservoir conditions of this 
facility are stagnation pressures of 2.41 to 10.0 MPa (350 to 
1450 psia) with a stagnation temperature of 1000 K (1800°/?) that 
produce perfect-gas freestream conditions of Mach 9.7 to 10.0 at unit 
Reynolds numbers of 0.5 x 10 6 to 2.0 x 10 6 /ft (1.6 x 10 6 1/m to 
6.6 x 10 6 /m). 

CEV compression pad model heating data were obtained in the 
Langley Research Center 31 Inch Mach 10 Air Tunnel during Test 



Fig. 7 Langley Research Center 31 Inch Mach 10 Air Tunnel. 


456. Test parameters were: angles of attack (a) of 20, 24, and 28 deg 
and freestream unit Reynolds numbers (Re 0 c ) of 0.5 x 10 6 to 
1.2 x 10 6 /ft. Nominal conditions for each of the test Reynolds 
numbers are given in Table 1; run-to-run variations in conditions 
were less than ±1%. For all test conditions, laminar flow was 
produced upstream of all cavities. The Mach 10 test data obtained 
during this study on the C2-D1 model supersede the C2-D1 data 
obtained during the previous [2] Test 446 (under the old IX model 
designation). This prior dataset is considered nonoptimal because the 
model support system had not been designed for that tunnel and 
positioning of the model to obtain an acceptable viewing angle and 
sufficient illumination was difficult. 

2. NASA Langley Research Center 20 Inch Mach 6 Air Tunnel 

The NASA Langley Research Center 20 Inch Mach 6 Air Tunnel 
(Fig. 8) is a blowdown facility in which heated, dried, and filtered air 
is used as the test gas; a detailed description of this facility can be 
found in [5]. The tunnel has a two-dimensional contoured nozzle that 
opens into a 20.5 x 20.0 in. (0.52 x 0.508 m) test section. The 
tunnel is equipped with a bottom-mounted injection system that can 
transfer a model from the sheltered model box to the tunnel centerline 
in less than 0.5 s. Run times of up to 15 min are possible in this 
facility, although for the current aeroheating study, run times of only 
a few seconds were required. The nominal reservoir conditions of 
this facility produce perfect-gas freestream flows with Mach 
numbers between 5.8 and 6.1 and unit Reynolds numbers of 0.5 x 
10 6 /ft to 7.3 x 10 6 /ft (1.64 x 10 6 to 23.3 x 10 6 /m). 

CEV compression pad model heating data were obtained in the 
Langley Research Center 20 Inch Mach 6 Air Tunnel during Test 
6954. Test parameters were: angle of attack (a) of 24 deg and 
freestream unit Reynolds numbers ( Re <*,) of 1.0 x 10 6 to 
6.9 x 10 6 /ft. Nominal conditions for each of the test Reynolds 
numbers are given in Table 2; run-to-run variations in conditions 
were less than ±1%. These conditions produced laminar flow 
upstream of the cavities for Reynolds numbers up to 5 x 10 6 /ft. At 
the Re ^ = 6.9 x 10 6 /ft condition, transitional flow was produced 
upstream from the three leeside cavities (pads nos. 6, 1, and 2). The 
C2-D1 and C3-D1 M 6 data from the previous Test 6946 [2] 
compression pad study (then with IX and 2X model designations) 
are also incorporated into the current analysis. 


Table 1 Nominal conditions for Langley Research Center 31 Inch Mach 10 Air Tunnel Test 456 


a, deg 

Re M , 1/ft 

Re oo, 1/m 


Poo, Pa 

Too* K 

Poo, kg/m 3 

Loo, m/s 

L 0 -// 30 ok, J/kg 

*FR,kg/m 2 /s 

<7fr, W/cm 2 

20, 24, 28 

5.19E + 05 

1.70E + 06 

9.68 

68.2 

52.4 

4.54E - 03 

1404 

7.38E + 05 

6.41E-02 

4.73 

20, 24, 28 

7.01E + 05 

2.30E + 06 

9.73 

93.9 

53.3 

6.15E-03 

1424 

7.66E + 05 

7.58E - 02 

5.80 

20, 24, 28 

1.02E + 06 

3.34E + 06 

9.81 

130.1 

51.9 

8.74E - 03 

1415 

7.53E + 05 

8.97E - 02 

6.75 

20, 24, 28 

1.22E + 06 

4.00E + 06 

9.84 

159.0 

52.7 

1.05E-02 

1431 

7.77E + 05 

9.98E - 02 

7.74 
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Fig. 8 Langley Research Center 20 Inch Mach 6 Air Tunnel. 


C. Data Acquisition, Reduction, Mapping, and Uncertainties 

1. Data Acquisition and Reduction 

Heating levels over the model were measured using the two-color, 
relative-intensity, global thermographic phosphor method [6,7]. In 
this method, heat- transfer coefficients are determined by assuming a 
step-function in heat transfer beginning at injection of the model into 
the tunnel, which corresponds to a parabolic temperature-time 
history. The model is illuminated by ultraviolet light sources that 
produce temperature-dependent fluorescence of the phosphor 
coating and images of the model are taken in the tunnel before and 
during a run using a three-color, charge-coupled device camera. The 
imaging for hypersonic experimental aerothermodynamic testing 
(IHEAT) code uses calibrations to convert the intensity data from 
each image pixel to temperatures and then performs the heat-transfer 
computations. Heat- transfer distributions from IHEAT are expressed 
in terms of the ratio h/h F R , where h FR is the heat- transfer coefficient 
resulting from a Fay-Riddell computation [8] for the heat transfer to 
a sphere with the radius of curvature of the CEV heat shield 
(8.4 in./0.21336 m) at a wall temperature of 300 K. 

2. Data Mapping 

The image data obtained from IHEAT were corrected for optical 
perspective effects and mapped to a three-dimensional surface model 
of the CEV geometry with compression pads. To accomplish this 
mapping, perspective transformations are first performed on the 
three-dimensional surface geometry until its two-dimensional 
projection matches that of the two-dimensional image data. The 
image data are then assigned transformed (. x , y, z) coordinates based 
on interpolation between the image and surface geometry, and then 
the transformation is inverted to obtain a corrected three-dimensional 
heating distribution. One refinement to this process has been made 
since its original use in the earlier compression pad test [2] . Instead of 
interpolating an (x, y, z) coordinate onto each image pixel, the image 
pixel heating data are instead interpolated onto each (. x , y, z) grid 
point of the surface geometry. While this method introduces some 
additional error into the interpolation of the image data, which is of 
much lower resolution than the surface geometry, the final result does 
not produce jagged edge at the bounds of the image domain and also 
retains the exact representation of the compression pad geometry. 



Fig. 9 Sample three-dimensional mapped image data, using inter- 
polation of image pixels to surface geometry. 


The differences between these two methods can be seen through 
comparisons of Figs. 9 and 10. 

After mapping the image data to the three-dimensional surface 
geometry, further information can be obtained by extracting heating 
distributions along specific line-cut stations. Data were extracted 
along the pitch- symmetry plane (i.e., the centerline) and along line 
cuts through the center of each pad that were roughly aligned with the 
near-surface streamlines at a = 24- deg. These compression pad 
data-cut locations are shown in Fig. 11. 

3. Data Uncertainties 

For testing of a blunt-body configuration such as the CEV, an 
uncertainty estimate for the measurement technique of approx- 
imately ±10% can be made based on [7]. Additional uncertainties 
are introduced through the mapping and processing of the data. 
Uncertainties also are introduced through the variation of wind- 
tunnel flow conditions. Since the repeatability of flow conditions is 
very good, the resulting uncertainty in heating measurements was 
estimated to be only on the order of ±5%. The other major source of 
uncertainty that must be considered is the mapping of pixel values 
from the image data to physical coordinates. This uncertainty is 
highly problem-dependent: on the overall acreage of a CEV model 
this uncertainty is estimated to be on the order of ±10%. However, in 
geometrically-complex regions, such as within/around the com- 
pression pads or at the heat-shield shoulder, where the scale of such a 
feature approaches the resolution of the camera optics, this 
uncertainty could be much higher. For instance, the localized heating 
spike observed downstream of a recessed pad’s cavity or on the 
leading edge of a protruding pad may be resolved to only a few image 
pixels and the resulting geometric uncertainty could be on the order 
of ±25 to ±50%. In the current analysis, the quoted uncertainty for 
the smooth OML was determined from a root-mean- square 
summation of that from the acquisition uncertainty (±10%), flow 
conditions (±5%), and the acreage estimate for image mapping 


Table 2 Nominal conditions for Langley Research Center 20 Inch Mach 6 Air Tunnel Tests 6946 and 6954 


a, deg 

Re 0 o, 1/ft 

Re^, 1/m 

Moo 

Poo- Pa 

T oo-K 

Poo- kg/m 3 

U oo , m/ s 

Hq—H^qo k, J/kg 

h FR , kg/m 2 /s 

q FR , W/cm 2 

24 

1.07E + 06 

3.5 IE + 06 

5.89 

295.3 

61.5 

1.67E-02 

926 

1.90E + 05 

7.48E - 02 

1.42 

24 

2.08E + 06 

6.83E + 06 

5.96 

576.3 

62.2 

3.23E — 02 

942 

2.05E + 05 

1.06E-01 

2.18 

24 

2.94E + 06 

9.63E + 06 

5.99 

803.0 

61.8 

4.53E - 02 

943 

2.06E + 05 

1.26E-01 

2.58 

24 

4.05E + 06 

1.33E + 07 

6.01 

1101.0 

61.7 

6.23E - 02 

844 

2.07E + 05 

1.48E-01 

3.05 

24 

5.04E + 06 

1.65E + 07 

6.02 

1415.1 

63.1 

7.83E-02 

957 

2.21E + 05 

1.68E — 01 

3.71 

24 

6.91E + 06 

2.27E + 07 

6.04 

1928.7 

63.0 

1.07E-01 

958 

2.21E + 05 

1.97E-01 

4.36 

24 

7.46E + 06 

2.45E + 07 

6.04 

2058.0 

62.5 

1.15E-01 

954 

2.17E + 05 

2.03E - 1 

4.40 
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Fig. 10 Sample three-dimensional mapped image data, using 
interpolation of surface geometry to image pixels. 


(±10%), which resulted in an overall experimental uncertainty of 
±15%. For the heating data in and immediately around the 
compression pads, the nominal mapping error was set to ±25% 
resulting in the overall pad-region uncertainty increasing to ±27%. 


D. Computational Method 

Flowfield predictions were performed using the Langley aero- 
thermodynamic upwind relaxation algorithm (LAURA) code [9,10] 
to obtain boundary-layer properties (e.g., Re e , M e , 8) for use in 
correlating the experimental data. LAURA is a three-dimensional, 
finite volume solver that includes perfect-gas, equilibrium, and 
nonequilibrium chemistry models. In this study, the perfect-gas air 
model was used for the wind-tunnel predictions. Freestream con- 
ditions in the Langley Research Center wind tunnels did not vary 
significantly from run to run, so the nominal conditions in Tables 1 
and 2 were used, with the wall temperature set to a constant 300 K. 
The use of a constant wall temperature was acceptable because the 
heat-transfer coefficient varies only slightly over the range of wall 
temperatures produced in these facilities. Solutions were performed 



at all wind-tunnel test conditions on a smooth-OML (no compression 
pad cavities) grid. The smooth-OML cases were performed using the 
laminar, thin-layer Navier-Stokes formulation. This methodology 
has been shown to produce good comparisons with experimental data 
in previous CEV studies, e.g., [11-13]. For purposes of data corre- 
lation and comparisons, heating rates and boundary-layer parameters 
such as Re e were taken from the OML solutions at locations 
corresponding to the centers of each pad. 

III. Results and Analysis 

A. Global Effects on Heating 

Sample three-dimensional mapped, global heating images 
obtained from the Mach 10 and Mach 6 tunnel tests are shown in 
Fig. 1 2 in terms of h /h FR . These data are for the C 1 -D 1 model and the 
images are arranged in order of increasing freestream Reynolds 
numbers. Data similar in character but of different heating levels 
were obtained for the other model configurations, although all 
models were not tested at all Mach 6 and Mach 10 Reynolds 
numbers. 

In these images it can be seen that the flow approaching the cavities 
was laminar at all test conditions except at the highest Mach 6 test 
Reynolds number ( Re ^ = 6.9 x 10 6 /ft) where transitional flow 
occurred upstream of the three leeside cavities. In all cases, the 
compression pad cavities perturbed the incoming flow leading to 
increased heating downstream of the cavities. For the three lowest 
Reynolds numbers (all of which are Mach 10 cases), there appeared 
to be little change in the magnitude of augmented heating down- 
stream of the pads, which indicated that the flow remained laminar; 
for these three cases, the heating on the cavity floor was 
also noticeably lower than that upstream of the pad. For Re ^ > 
1.22 x 10 6 /ft, the three leeside pads (nos. 1, 2, and 6 at 0 deg and 
±60- deg) experienced increasing downstream heating with 
Reynolds numbers which indicated progression of the boundary- 
layer state at the pad from transitional to turbulent flow. For the no. 3 
and 5 pads (at ±120- deg) similar behavior was observed but with 
transition onset appearing first at the Mach 6, Re ^ = 2.08 x 10 6 /ft 
condition. Local Reynolds numbers at the no. 4 pad (near the 
stagnation region) remained very low and transition onset was not 
observed until the highest Mach 6 Reynolds number condition of 
Re OQ = 6.9l x 10 6 /ft. 

B. Detailed Cavity-Heating Distributions 

Heating distributions extracted from the image data in Fig. 12 
along the line cuts through each pad are shown in Fig. 13. For each 
pad, the data from all test Reynolds numbers are plotted vs z/R for 
pad nos. 1,2,4, and 6 or, to better reflect the orientation of the line cut, 
vs y/R for pads nos. 3 and 5. The upstream h/h FR distributions were 
nearly constant for all Reynolds numbers except at Re 00 = 
6.91 x 10 6 /ft, which provided evidence that the flow incoming to 
the pads was laminar for all conditions except at this highest 
Reynolds number where it was transitional upstream of pads nos. 1, 
2, and 6. Within the compression pad cavities, the h/h FR distri- 
butions also remained relatively constant at the three lowest 
Reynolds numbers, again indicating laminar flow. At Re ^ > 
1.22 x 10 6 /ft (depending on the pad location), the h/h FR distri- 
butions began to increase with Reynolds number as is char- 
acteristic of transitional or turbulent flow. For all cases, the highest 
heating levels were produced around the downstream wall and lip of 
each cavity, which essentially represents a forward facing step. 

C. Cavity-Heating Augmentation 

In the design of the CEV flight vehicle, the heating augmentation 
above smooth-OML levels in and around the compression pad 
cavities is of interest; this quantity will be referred to as h/h SM00TH . 
To determine this ratio, a smooth-OML reference heating value, 
(h/h FR ) SM ooth> was first determined for each pad at each test 
condition by averaging the measured heating levels over a small 
region upstream of the respective pad. The augmentation factor 
h/h s mooth was then computed by dividing the measured h/h FR 
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M 0 o = 10, Re ao = 0.52x1 0 6 /ft 


Moo = 10, Re«> = 1.22xl0 6 /ft 





Moo =10, Re x = 0.70x1 0 6 /ft 


Moo = 6, tooo = 2.08x1 0 6 /ft 


Moo = 6, Re«,= 5.04x1 0 6 /ft 




Moo =10, Re x = 1.02xl0 6 /ft 


Moo = 6, Re oo = 2.94xl0 6 /ft 




Moo = 6, Re* — 6.91x 10 6 /ft 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Fig. 12 Global heating distributions measured on Cl-Dl model, a = 24- deg. 


values around each pad by the reference (h/h FR ) smooth f° r that pad. 
The heating data from Fig. 12 is shown recast in terms of h/ ^smooth 
in Fig. 14. 

The augmentation factors were — 1 around the pads, or <1 within 
the pads, for the lowest Reynolds numbers at which conditions the 
flow in and downstream of the pads remained laminar. For Re ^ > 
1.22 x 10 6 /ft, augmentation factors increased with Reynolds 
numbers to maximum values around 3 as the flow progressed from 
transitional to turbulent in and around the pads. Close examination of 
the augmentation factors around pads nos. 1, 2, and 6 (0 and 
±60- deg) shows that the augmentation factors actually decreased at 
the highest test Reynolds number of Re ^ = 6.91 x 10 6 /ft. The 
decrease was due to the definition of the augmentation factor in this 
analysis, which is the measured pad heating relative to the measured 
upstream heating. At this highest Reynolds number, the flow 
upstream of these pads was already transitional and thus the 
measured reference values were higher than if the flow had remained 
laminar. So, although the actual dimensional heating was still the 
highest at this condition, the relative augmentation values were 
actually slightly lower than those of the next lowest Reynolds 
number conditions. 

D. Cavity Transition and Average Heating 
Augmentation Correlations 

1. Range of Cavity Average Heating Augmentation Data 

Another parameter of interest to design of the CEV flight vehicle is 
the average heating augmentation over the surface top of each pad 


(i.e., the cavity floor), which is denoted as {h/h smooth) avg- These 
values were determined for all pads and test conditions and the 
resulting values are plotted in terms of the boundary-layer momen- 
tum thickness Reynolds number, Re e , in Fig. 15. These average 
augmentation factors varied from —0.5 to —3.0, however, at first 
glance there seemed to be no correlation with Re e . Although not 
shown, a similar lack of correlation was found when the data were 
plotted against Re e /M e . The maximum augmentation factors both 
on and downstream of the pad are also quantities of interest. 
Maximum values varied between 1 and 4, but the analysis of those 
data has not yet been completed. 

More insight into the physics of this problem can be obtained when 
the data are color-coded with respect to the state of the boundary layer 
at and upstream of the compression pad cavities (Fig. 16), as 
determined by examination of the data plotted and displayed as in 
Figs. 12-14. For conditions where the incoming flow was laminar 
and remained laminar downstream of the cavity, the average 
augmentation (/*//zsmooth)avg decreased with Re e to a minimum 
value of {h/h smooth) avg ~ 0.4. For cases where the cavity caused 
transition, average augmentation factors increased to a maximum of 
(^Asmooth)avg ~ 2.5. And finally, for cases where the incoming 
flow was transitional or turbulent, the average augmentation factor 
decreased (as defined relative to measured upstream levels) from the 
2.5-maximum value. Of particular interest in this plot is that the 
average augmentation factor at the pad approaches 1 at the highest 
Reynolds numbers, which indicates that for fully-turbulent incoming 
flow conditions, the average pad heating is no higher (approximately) 
than the incoming turbulent level. 
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0 n — 1 — « — I ■ 1 ■ ■ ■ ■ I ■ ; 1 — L 1 

0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00 

z/R 


-0.75 


-0.70 


-0.65 


-0.60 


Mach 10, Re -0.52*1 0*/H 
Mach 10, Re.«0.70*10‘/tt 

Mach 10, Re -1.02*1 OVft 

— • — Mach 10, Re -1.22x1 o'/lt 

• Mach 6, Re. -2.08x1 oVft 

■ Mach 6, Re. -2.94x1 oVft 

• Mach 6, Re. -4.05x1 oVft 

• Mach 6, Re -5.04x1 0*/tt 

■ Mach «, Re.-6.91 -10*711 


4.o r 


3.5 - 


3.0 - 
2.5 - 


oc 


^ 2.0 - 

.c 


1.5 - 


1.0 - 


0.5 - 
°-8.60 


z/R 



_J I I I I I I I 

0.55 0.50 0.45 0.40 0.35 0.30 0.25 0.20 



z/R 



Fig. 13 Heating distributions through Cl-Dl model pad centers, a = 24- deg. 


These data clearly indicate that the cavity effects are a function of 
the boundary-layer state. Thus, in order to develop a correlation for 
the heating augmentation factors, a correlation must first be 
developed that can be used to define the onset of transition at a 
compression pad cavity. 

2 . Crew Exploration Vehicle Closed-Cavity Transition Correlation 

A first attempt at correlating cavity-induced transition on the CEV 
was performed by applying transition criteria developed for cavities 


on the Mars Science Laboratory (MSL) heat shield [14-16]. 
Although the carrier-vehicle attachment points were eventually 
eliminated from the MSL heat shield by moving them to the aftbody 
of the vehicle, an extensive database was developed through wind- 
tunnel testing and analysis that also included an earlier similar study 
for the Genesis mission [17]. This database was used to develop the 
correlations shown in Fig. 17 for cavity size and flow conditions that 
bound fully-laminar and fully-turbulent flow produced at a cavity. 

Transition-onset values were determined from the current CEV 
data set, as well as from the data in [3], which provided smaller cavity 
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Fig. 14 Heating augmentation factors around Cl-Dl model pads, a = 24- deg. 


aspect ratio data {D/H = 5 to 10) as well as data {D/H = 1 1 to 18) 
that overlap the current data set {D/H = 1 1 to 45). These values are 
compared with the MSL laminar/transitional boundary in Fig. 18, 
which is comprised almost entirely of data points from cavities with 
aspect ratios in the D /H = 1 to 4 range. Unfortunately, the only CEV 


data that come close to matching the MSL correlation are those for 
the D/H = 5 to 6 data; the remaining CEV data clearly diverge 
farther from the MSL curve as a function of the cavity aspect ratio. 

The reason that the MSL correlation does not work for the CEV 
data is that because of the difference in sizes between the two sets of 
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Fig. 15 Average heating augmentation factors for all pads, all test Fig. 16 Average heating augmentation factors classified by boundary- 
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Fig. 17 Laminar and turbulent boundaries from MSL transition study 
[16]. 


cavities, the CEV and MSL cavities flows represent different flow 
regimes. As shown in Fig. 19, a delineation is typically made (e.g., 
[18,19]) between small aspect ratio {D/H < 10) and large aspect 
ratio {D/H >15) cavities. An open cavity (e.g. the MSL data) 
produces separated, recirculating flow with the cavity, whereas the 
flow attaches to the cavity floor in a closed cavity (e.g. the majority of 
the CEV data). 

Following the definitions presented in Fig. 19, the C0-D1 and 
Cl-Dl models {D/H = 45.6 and 22.8, respectively) can clearly be 
expected to produce closed-cavity, attached flow, whereas the C3-D1 
model {D/H = 5.7) should produce open-cavity, separated flow. 
The C1-D2 and C2-D1 models (both with D/H = 1 1 .4) fall into the 
intermediate cavity regime, but will provisionally be treated as closed 
cavities. Additional CEV data from [3] span the range from open to 
closed-cavity flow. 

With respect to the actual CEV flight vehicle, current designs call 
for very large aspect ratio (wide and shallow) cavity design. Thus, a 



Fig. 18 CEV transition-onset data compared with MSL laminar 
boundary criterion. 



H/6 

Fig. 20 CEV transition criteria for closed-cavity, attached flow. 


transition correlation relevant to closed-cavity, attached flow is 
required. Because the MSL methodology (in which transition was 
shown to correlate in terms of Re e and the ratio of cavity diameter to 
boundary-layer height, D/8) was shown to be invalid for CEV, a 
different functional relationship was required. Since closed-cavity 
flows produce attached flow on the cavity floor and then the cavity 
end wall provides the flowfield disturbance that causes transition, an 
analogy to transition correlations for discrete boundary-layer trips 
can be made. For such correlations (e.g. [20]), boundary-layer 
transition is treated as a function of Re e /M e (rather than just Re e ) and 
the ratio of trip height (or in this case cavity depth) to boundary -layer 
thickness, H/8. 

Such a correlation is presented in Fig. 20 based on all CEV data for 
11.5 < D/H < 45.6. A good fit was produced with a correlation 
coefficient of 0.86. The divergence of the open-cavity, separated flow 
CEV, MSL and Genesis data from this fit can also be seen. If the 
intermediate cavity data {D/H = 11.5, which are not relevant to 
current CEV designs) are excluded, the correlation coefficient can be 
increased to 0.97 as shown in Fig. 21 . So, depending on the range and 
accuracy required, either of the following relations made be used to 
predict the onset of transition due to compression pad cavities for 
closed-cavity, attached flows: 



Having developed a boundary-layer transition criteria, separate 
correlations for laminar-in/laminar-out and laminar-in/transitional- 
or-turbulent-out average heating augmentation factors now can be 
developed. However, a few provisions on the applicability of all these 
correlation must first be provided. At the test conditions of this study, 
the edge Mach numbers on the blunt forebody of the CEV models 
ranged from M e = 0.1 to 1.1; thus there is little overlap between 


Open-cavity, separated flow 
D/H <10 (approx.) 





Intermediate cavity Closed-cavity, attached flow 

1 0 < D/H <15 (approx.) 1 5 < D/H (approx.) 



Fig. 19 Illustration of open and closed-cavity flow regimes. 
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Fig. 21 CEV closed-cavity, attached-flow transition criteria refined to 
large ( D/H ). 


these data and data on flat plates or lifting vehicles such as the shuttle, 
X-33, or X-38 with edge Mach numbers ranging from — 1 to —4. 
Also, these data were obtained in perfect-gas wind tunnel in 
which the parameters 8/0 and T w /T e , which are sometimes em- 
ployed in transition correlations, varied only slightly. Approxi- 
mate relations for these parameters as a function of Re e are provided 
in Figs. 22 and 23. 


3. Crew Exploration Vehicle Average Cavity-Heating 
Augmentation Correlations 

The pad-averaged heating augmentation factor data that corre- 
spond to laminar-in/laminar-out conditions, as identified using the 
transition criteria of Eq. (1) or Eq. (2), are plotted in Fig. 24 vs Re d . 
These data were sorted by the cavity aspect ratio and showed a clear 
trend of increasing augmentation with increasing aspect ratio; or, 
looking at it from the other direction, the cavity floor became cooler 
as the cavity became deeper. Based on this trend, a correlation 
variable was defined by dividing Re e by the cavity aspect ratio, as 
shown in Fig. 25. Plotted in this matter, the data showed a clear 
“^-curve” shape with augmentation factors that varied from a 
maximum of —1 to a minimum of —0.3. This behavior is similar to 
that of a hyperbolic tangent function, and so a correlation equation of 
that form was developed for these data. This laminar-in/laminar- 
out closed-cavity attached-flow augmentation factor correlation is 
given by 



Fig. 22 Approximate relation for variation of 8/0 as a function of Re 0 
on the CEV models. 



Fig. 23 Approximate relation for variation of T w /T e as a function of 
Re e on the CEV models. 

r= (l — 'I x,anl P' 718 < X - M2) l 

V Ml) \Ml) L M 2 

Ml =4.0152, M2 = 2.3589 

where X = Re e x ( — \ and Y = ( — ^ (3) 

V^SMOOTH/AVG 

This function has a correlation coefficient of 0.787 and fits the 
majority of the C0-D1 and Cl-Dl model cavity data (D/H of 45.6 
and 22.8, respectively) to within ±15% limits. The major 
divergences from this curve fit were from the C1-D2 and C2-D1 
model data; these two models both had aspect ratios of D /H = 11.5, 
which are between the open and closed-cavity limits. There was also 
fair amount of scatter at the very low X values where the augmen- 
tation should ideally approach 1. These data were mostly from pad 
no. 4, which is in the stagnation region. These data were harder to 
optically acquire than elsewhere on the model because the stagnation 
region surface was more inclined away from the camera line of sight 
due to curvature of the model. Also, the boundary-layer variables 
extracted from the numerical computations are thought to be less well 
defined since the flow velocity approaches zero at the stagnation 
point. 

The laminar-in/transitional-or-turbulent-out heating augmenta- 
tion factor data was found to require a different correlation function 



Fig. 24 Laminar, closed-cavity, attached-flow augmentation factors vs 
Re e . 
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Fig. 25 Correlation function for laminar, closed-cavity, attached-flow Fig. 27 Transitional/turbulent, closed-cavity, attached-flow augmen- 
augmentation factors. tation factors vs Re e /M e . 


than the laminar-in/laminar-out data. A plot of these data is shown in 
Fig. 26 in terms of Re e with the values color-coded by pad 
aspect ratio. Three trends were evident from this plot: first, that 
the augmentation factors increased with Re e \ second, that the aug- 
mentation factors increased with pad depth; and finally, that the data 
from each pad on a given model fell into separate regions. 

To collapse these data, they were first plotted in terms of Re e /M e , 
as shown in Fig. 27; however considerable spread was still evident. 
Through a trial- and-error process, it was determined that plotting the 
data in terms of ( Re e /M e ) x (1 /M e ) l/2 provided a much better 
collapse as shown in Fig. 28, although a dependence on aspect ratio 
was still evident. Additional trial-and-error fitting led to the inclusion 
of the factor (D/H)~ 024 and a final correlation function, as shown in 
Fig. 29, of the form 


Y = Ml x X M1 , 


Ml = 0.003574, 


M2 = 1.2443 


wi, " ex =(t ! )(ii) l, ’(s) 054 and 

y= — ) 

V^SMOOTH/ avg 


( 4 ) 


As with the laminar augmentation correlation, this function fits the 
majority of the high cavity aspect ratio data from the C0-D1 and 


Cl-Dl models to within ±15% limits. Greater differences were seen 
in the C2-D1 and C3-D1 model data, which have cavity aspect ratios 
of 11.5 and 5.7. 

4. Crew Exploration Vehicle Open-Cavity Correlations 

At this point, no attempt has been made to develop transition or 
heating augmentation correlations that encompass both open and 
closed-cavity data since open-cavity flows are not likely to be 
experienced by the CEV vehicle. However, such correlations for 
blunt-body, subsonic boundary-layer flow may be possible (if more 
complex) and will be explored in the future. One such unified open 
and closed-cavity correlation for supersonic boundary-layer-edge 
flow has recently been developed for Shuttle Orbiter damage 
assessment [21,22]. 

5. Turbulent and Tripped Data Correlations 

The correlations developed herein do not take into account data 
from the wind-tunnel tests in which the flow entering the com- 
pression pad cavities was already turbulent, whether from natural 
boundary-layer transition or due to discrete trips placed upstream of 
the cavities. While these data will be examined, it was shown in 
Fig. 16 that the augmentation factors (as defined relative to the actual, 
measured upstream heating) for such cases were actually lower than 




Fig. 26 Transitional/turbulent, closed-cavity, attached-flow augmen- Fig. 28 Transitional/turbulent, closed-cavity, attached-flow augmen- 
tation factors vs Re e . tation factors vs ( Re e /M e ) (1/M e ) 1/2 . 
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X = (Re/M e )(1/M e ) 1/2 (D/H) 0 - 24 


Fig. 29 Correlation function for transitional/turbulent, closed-cavity, 
attached-flow augmentation factors. 

for cases with laminar incoming flow followed by transition 
produced by the cavity. 

IV. Conclusions 

The effects of recessed compression pads on the aeroheating 
environment of the CEV heat shield have been studied through 
testing of a CEV models with various cavity dimensions at Mach 6 
and 10. Global heating data were obtained and mapped to three- 
dimensional surface geometries and detailed heating distributions 
were extracted from these maps along line cuts through the center of 
each pad. Augmentation factors for the cavity effects on heating were 
determined by taking the ratio of measured heating at and around the 
pad to that measured on the smooth OML upstream of the pad. 

The test conditions provided laminar incoming flow to all pads 
except for the three leeside pads at the highest test Reynolds number 
at Mach 6. Depending on pad geometry and local conditions, the 
boundary-layer state at and downstream of the pads varied from 
laminar to transitional to turbulent. Compression pad cavity floor- 
averaged augmentation factors varied from a minimum of approxi- 
mately 0.3 at laminar-in/laminar-out conditions to a maximum of 
approximately 2.5 at laminar-in/turbulent-out conditions. 

A cavity transition correlation was developed in terms of the 
boundary-layer parameter Re e /M e and the ratio of cavity depth to 
boundary-layer height, H/8. This correlation is valid for closed- 
cavity, attached flow with cavity aspect ratios of D/H >10, which 
are applicable values to the design of the CEV vehicle. It was also 
shown that open-cavity, separated flow ( D/H < 10) transition data 
do not fit this correlation; however such values are not within the 
CEV vehicle design space. 

This transition correlation was used to classify the measured 
augmentation factors as either laminar or transitional/turbulent and 
separate correlation functions were developed for each data set. A 
majority of the closed-cavity heating augmentation factor data fell 
with ±15% of these correlation fits, while larger differences were 
evident in the open-cavity data. 
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